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Abstract Bats have great potential as reservoirs for 
emerging viruses such as severe acute respiratory syn- 
drome-coronavirus. In this study, bat coronaviruses 
(BtCoVs) were detected by RT-PCR from intestinal and 
fecal specimens of Miniopterus fuliginosus breeding colo- 
nies in Wakayama Prefecture caves, where we previously 
identified bat betaherpesvirus 2. Two primer sets were used 
for the detection of BtCoV: one was for the RNA-depen- 
dent RNA polymerase (RdRp) region and the other was for 
the spike (S) protein region. Eleven and 73% of intestinal 
and fecal specimens, respectively, were positive for RdRp 
region, and 2 and 40% of those were positive for S protein 
region. Sequencing and phylogenetic analysis showed that 
the detected BtCoV belonged to the group 1 (alpha) cor- 
onaviruses. These data suggest that BtCoV is endemic in 
M. fuliginosus in Japan. 
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Introduction 


An outbreak of severe acute respiratory syndrome (SARS) 
that was caused by a newly identified SARS coronavirus 
(SARS-CoV) occurred in China, Singapore, Vietnam, and 
other countries from 2002 to 2003; 8,098 patients were 
reported, 774 of whom died [1]. At first, SARS-CoV was 
considered to be derived from civets and raccoon dogs [2]. 
However, recent studies have suggested that SARS-CoV is 
a recombinant of bat-derived coronaviruses [3, 4]. Bat 
coronaviruses (BtCoVs) have now been identified in China 
(3, 4] and other countries such as USA [5], Germany [6], 
Kenya [7], and the Philippines [8]; however, there has been 
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no report of BtCoV detection in Japan. Bats have also great 
potential as reservoirs for emerging viruses such as Ebola 
and Nipah [9, 10]. We have recently identified novel 
viruses from bats, such as bat betaherpesvirus 2 [11] and 
bat adenovirus | [12], in Japan. To confirm the role of bats 
as host species for viruses, we attempted to detect BtCoV 
from Miniopterus fuliginosus using intestinal and fecal 
specimens from specific caves located in the Wakayama 
prefecture, Japan. 


Materials and methods 
Bat samples 


All bats and fecal specimens were collected in May 2009 and 
July 2010 from a headrace tunnel and a cave in the Wakay- 
ama Prefecture, Japan, with permission from government 
officials. Those are well-known habitats of M. fuliginosus; 
the headrace tunnel is a roosting place of male and non- 
breeding female bats and the cave is a breeding place, where 
we collected bats and identified a novel virus in a previous 
study [11, 13]. The following three methods were used to 
collect intestinal specimens: (1) Captured bats were eutha- 
nized using an overdose of anesthetic and small intestines 
were collected without washing their contents, (2) Forty-five 
bats were captured with a net and placed in a bag, after which 
fresh fecal samples were collected, placed in a solution of 
RNAlater (Ambion, Foster City, CA, USA), and separated 
into four pools. After fecal samples were collected, the bats 
were released safely, (3) Three traps were placed under a 
mass of bats that were hanging from the ceiling of each cave. 
The next day, fecal samples were collected, placed in 
RNAlater, and separated into 11 pools. 


Extraction of RNA and RT-PCR 


Total RNA was first extracted using TRIzol (Invitrogen, 
Carlsbad, CA, USA) and was re-extracted using the 
QIAampViral RNA Mini Kit (Qiagen, Hilden, Germany). 
The cDNA was synthesized using M-MLV_ Reverse 
Transcriptase (Invitrogen, San Diego, CA, USA) and ran- 
dom hexamers, and two-step reverse transcription PCR was 
performed using PrimeSTAR Max (Takara Bio, Shiga, 
Japan) with the default conditions of each master mix, 
except the annealing temperature, which was 48°C. Two 
primer sets were used for the detection of BtCoV; primer 
set for the first PCR: IN-6 primer, 5’-GGTTGGGAC 
TATCCTAAGTGTGA-3’ and IN-7 primer, 5’-CCAT 
CATCAGATAGAATCATCAT-3’ [14]; primer set for the 
second PCR: IN-6 primer and hemi-nested reverse primer, 
5'-ATCAGATAGAATCATCATAGAGA-3’ [15], were 
used for detection of the sequence in RNA-dependent RNA 
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polymerase (RdRp) region. These primer sets are known to 
be able to amplify genes of the coronavirus family 
including BtCoV [5, 8, 14, 16-18]. To amplify the spike 
(S) region sequence, consensus primers were constructed 
using the online consensus primer design software 
(CoCoMo; http://www.geneknot.jp/cocomo/) with seven 
sequences deposited in GenBank. The constructed primer 
sets for S protein region were as follows: No. 6, 5’- 
NSHRYKTATGTHTGYAAYGGHAA-3’ and No. 2 5’- 
DGAYTGBGAYTTDACACAYTCRTT-3’ were used for 
the first PCR; No. 1, 5'-HTGTGYBCAGYAYTAYAAY 
GGYAT-3’ and No. 2 were used for the second PCR. After 
amplification, amplicons were detected by 1.5% agarose 
gel electrophoresis. 


Phylogenetic analysis 


The phylogenetic analysis was performed as follows: 
obtained sequence data were carefully checked for peaks 
and 327 bp (RdRp) and 547 bp (S) of each amplicon, 
except primer sequences were selected to construct phy- 
logenetic trees of nucleotide sequences. Five different 
sequences were obtained for the RdRp region (accession 
numbers AB619638—AB619642) and two sequences were 
obtained for S region (accession numbers AB644273 and 
AB644274). The same region of other available coronav- 
iral sequences (BtCoV HKU7-1, DQ249226; BtCoV 
HKU8-1, DQ249228; BtCoV 1A, NC_010437; BtCoV 1B, 
NC_010436; BtCoV Fujian/773/2005, EF434379; HCoV 
229E ATCC, NC_002645; HCoVNL63, NC_005831; 
PEDV CV777, NC_003436; AIBV Beaudette GK, 
AJ311317; MHV JHM, NC_006852; SARSr-Rh-BatCoV 
HKU3-7, GQ153542; SARS Frankfurtl, AY291315; civet 
SARSr-CoV PC4-139, AY613949; and EU420139, BtCoV 
HKU8) were used to create the phylogenetic trees. The 
bootstrap test for the phylogenetic tree was constructed 
using MEGA4 software [19] with the neighbor-joining 
method with a pairwise deletion option. 


Results 


As shown in Table 1, PCR products at the expected sizes 
were obtained from 5/45 (11%) of intestinal and 11/15 
(73%) of fecal specimens with the RdRp primer set, and 
1/45 of intestinal and 6/15 (40%) of fecal specimens were 
positive with the S protein primer set. PCR products were 
gel purified and sequence analysis was performed with 
specific-sequence primer. The sequences of PCR products 
showed some similarities to the deposited BtCoV sequen- 
ces; these sequences were analyzed by a BLASTn search 
(http://blast.ddbj.nig.ac.jp/top-e.html) and matched to 
nucleotide sequences of BtCoVs. These suggest that the 
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Table 1 Detection of BtCoV from bat specimens 
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Region Specimens RT-PCR Sequence type Number Total 
RdRp Intestine 5/45 (11%) BtCoV M.ful./Japan/01/2009 4 5/45 (11%) 
BtCoV M.ful./Japan/02/2009 1 
Feces from 
Captured bat 2/4 (50%) BtCoV M.ful./Japan/01/2010 1 11/15 (73%) 
BtCoV M.ful./Japan/03/2010 1 
Trap 9/11 (82%) BtCoV M.ful./Japan/01/2010 7 
BtCoV M.ful./Japan/02/2010 2 
S Intestine 1/45 (2%) BtCoV M.ful./Japan/03/2009 1 1/45 (2%) 
Feces from 
Captured bat 1/4 (25%) BtCoV M.ful./Japan/04/2010 1 6/15 (40%) 
Trap 5/11 (45%) BtCoV M.ful./Japan/04/2010 4 


detected sequences were derived from newly identified 
BtCoV sequences. Because the PCR product was detected 
not only in the feces pool but also in the small intestines, 
these results strongly suggest that M. fuliginosus in Japan is 
definitely infected by BtCoV. 

Obtained sequence data were carefully checked for 
peaks and 327 bp for the RdRp region and 547 bp for the S 
protein region of each amplicon, except primer sequences, 
were selected to construct phylogenetic trees of nucleotide 
sequences. Five different sequences of BtCoVs were 
detected in RT-PCR products of the RdRp region and two 
sequences were detected from RT-PCR products of the S 
protein region, and they were named in accordance with 
the rules for nomenclature of influenza virus, that is, 
BtCoV_ host species/isolated place/isolation number/ 
isolated year. The viral genomic information was deposited 
in DDBJ/EMBL/GenBank (accession numbers AB619638-— 
AB619642 and AB644273-AB644274). The phylogenetic 
trees were constructed along with other available coro- 
naviral sequences deposited in GenBank using the MEGA4 
software [19] (Fig. la: RdRp region, b: S region). The 
sequences detected from M. fuliginosus in Japan were 
classified into group | (alpha) coronavirus of BtCoVs and 
were quite different from SARS-related BtCoVs. In the 
RdRp region, BtCoVs in Japan are closely related to 
BtCoV HKU7-1, which was isolated from M. magnater in 
Hong Kong in 2006 [18], with 95% of nucleotides identi- 
fied. BtCoV HKU8-1, which was isolated from M. pusillus 
in Hong Kong in 2006 [18], was the second related strain, 
but fewer nucleotides were identified (80%) than in HKU7-1. 
In the S protein region, M.ful./Japan/03/2009 showed 92% 
identity relative to BtCoV Fujian/773/2005 (GenBank 
EF434379), which was isolated in the Fujian province in 
China in 2005 [20]. In contrast, M.ful./Japan/04/2010 
showed about 70% similarity with the HKU8 strain and 
other group 1 BtCoVs (1a and 1b), which were isolated in 
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Hong Kong from 2004 to 2005 [21]. Both of two isolated 
Japanese sequences belonged to the group 1 (alpha) coro- 
naviruses, but the identity between them was only 64%. 
Therefore, this suggests that different BtCoV strains were 
present in the same cave. We could not compare nucleotide 
sequences of the S gene between Japanese BtCoV and 
HKU 7-1 because the sequence for HKU 7-1 was not 
available in GenBank. 


Discussion 


In this study, 45 intestinal specimens and 15 pools of fecal 
specimens were tested for BtCoV with specific primer sets 
for RdRp and S regions. The number of PCR-positive pools 
was high: 5/45 (11%) of intestinal and 11/15 (73%) fecal 
specimens were positive with the primer set for the RdRp 
region and 1/45 (2%) of intestinal and 6/15 (40%) fecal 
specimens were positive with the primer set for the S 
protein region. The specimens were collected from two 
places; one was a habitat of M. fuligunosus, which lives in 
the Kinki region of Japan, and more than 20,000 female 
bats gather in the cave during the breeding season [13]. 
This suggests that BtCoV is endemic at a high frequency in 
M. fuliginosus in the Kinki region of Japan. In addition, 
BtCoVs in Japan belong to the group 1 (alpha) coronavi- 
ruses, and phylogenetic analysis of the RdRp and S protein 
regions indicated that it is quite different from SARS- 
related CoV. In particular, S protein of coronavirus is a 
major determinant of receptor binding and virus—cell 
membrane fusion [22]. Therefore, the differences in S 
protein region sequence imply the possibility that Japanese 
BtCoVs have different infectivity to host species compared 
to other SARS-related BtCoVs. Moreover, the S protein 
sequences of M.ful./Japan/03/2009 and M.ful./Japan/04/ 
2010 showed 64% identity between them. This also 
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Fig. 1 The phylogenetic trees (a) 
based on the nucleotide 
sequence of a RdRp and b S 
regions. The phylogenetic trees 
were constructed with the 
MEGA4 software using the 
bootstrap test command with the 
neighbor-joining method. The 
accession numbers of used 
coronaviral sequences were 
described in “Materials and 
methods” section 


43 


Group 


BtCoV M.ful./Japan/01/2009 AB619638 
BtCoV M.ful./Japan/02/2010 AB619641 


BtCoV M.ful./Japan/02/2009 AB619639 
BtCoV M.ful./Japan/01/2010 AB619640 
BtCoV M.ful./Japan/03/2010 AB619641 
BtCoV HKU7-1 DQ249226 
BtCoV Fujian/773/2005 DQ648835 I 
BtCoV HKU8-1 DQ249228 

BtCoV 1A NC_010437 

BtCoV 1B NC_010436 

HCoV229E ATCC NC_002645 
HCoVNL63 NC_005831 

PEDV CV777 NC_003436 

AIBV Beaudette CK AJ311317 
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0.1 


suggests that different strains of BtCoVs are endemic in the 
same place, like several different strains of human 
coronaviruses. 

The nucleotide sequences of RdRp of BtCoV in Japan 
were similar (95%) to those of HKU7-1, which was iso- 
lated from M. magnater in Hong Kong, and the nucleotide 
sequences of S protein of Japanese BtCoVs also showed 
some similarity to other Hong Kong and Fujian strains. 
Miniopterus bats are a migratory species and one popula- 
tion has several habitats surrounding its central breeding 
cave [23]. The distance of bat migration is reported to be 
several 100 km [13, 24]. This implies the possibility that 
BtCoV-infected bats have been brought to Japan from 
Southeast China (or to there from Japan) through unknown 
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routes. In addition, Rodrigues and Palmeirim [24] have 
reported that female bats nearly always return to their natal 
colony to give birth, whereas male bats sometimes go to 
other colonies, suggesting that BtCoV is transmitted 
between Japan and Hong Kong or Fujian by male bats. 
Similar cases have been reported for rabies virus. Rabies 
virus infection is associated with the migratory routes of 
Nathusius’ pipistrelle (Pipistrellus nathusii) in France. It 
has also been reported that silver-haired bats (Lasionycteris 
noctivagans) migrate seasonally from Alaska, across Can- 
ada, to Texas, and rabies virus variants have been identified 
from several locations throughout the geographic range of 
these bats [25]. Therefore, BtCoV might be transmitted 
between Hong Kong and Japan as a result of bat migration. 
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